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We point out that in the irreducible natural SUSY spectrum, stops have comparable branching fractions to 
chargino-bottom and neutralino-top in the vast bulk of parameter space, provided only that both decay modes 
are kinematically accessible. The total stop pair branching fractions into ti+ MET can therefore be reduced 
to 0(50%), while bb + X branching fractions are typically much smaller, 0(10%), thus limiting the reach of 
traditional stop searches. We propose a new stop search targeting the asymmetric final state it* — > tx°bx~ + 
H.c., which can restore sensitivity to natural stops in the 7 and 8 TeV LHC runs. In addition we present a 
new variable, topness, which efficiently suppresses the dominant top backgrounds to semi-leptonic top partner 
searches. We demonstrate the utility of topness in both our asymmetric search channel and traditional tt* — > 
tt+MET searches and show that it matches or out-performs existing variables. 

PACS numbers: 



INTRODUCTION 

The recent discovery of what looks very much like a 
weakly-coupled Higgs boson has made the question of the 
electroweak hierarchy acute. Weak-scale SUSY has long been 
a favored contender to protect the Higgs vacuum expectation 
value, but the impressive agreement of LHC data with SM 
predictions has made it clear that most of the MSSM parti- 
cle content must lie beyond the kinematic reach of the 7 and 
8 TeV LHC if i?-parity is conserved. Only a small subset of 
the SM superpartners are immediately relevant for the natu- 
ralness of the Higgs potential, however |fl]. In particular, only 
Higgsinos, with masses mjj 5s 200 GeV, and third-generation 
squarks tr,,tR, br,, with masses nig < 400-500 GeV, need to 
be light to keep EWSB fully natural, though gluinos, winos, 
and binos must not be too far beyond current LHC energies 
12 • Searching for a light stop, or three light third-generation 
squarks, therefore becomes critical for understanding the ex- 
tent to which SUSY is relevant for stabilizing the electroweak 
scale. 

The purpose of the present work is twofold. First, we point 
out that, in a natural SUSY spectrum, the pure Higgsino nature 
of the light neutralinos and chargino can significantly weaken 
the reach in the traditional stop search channel, tt* — > ti+I^x, 
and we propose a novel search channel targeting the natural 
SUSY spectrum which can recover the lost sensitivity. Sec- 
ond, we introduce a new kinematic variable, topness, which 
efficiently suppresses the dominant backgrounds to top part- 
ner searches in semileptonic final states. We demonstrate its 
performance in both our novel search topology and in more 
traditional searches. 



STOP AND SBOTTOM BRANCHING RATIOS 

In the natural SUSY spectrum, the lightest states 
are Higgsinos and third generation squarks; the gaug- 
inos can be heavier, up to several hundreds of GeV 
Mass mixing between the higgsinos and the heavier 



EW gauginos splits the Higgsinos; for example, with 
(Mi , M 2 , n, tan /?) = (400 GeV, 800 GeV, 200 GeV, 20), 
the resulting Higgsino masses are (m x o,m x ±,m x o) = 
(192 GeV, 197 GeV, 204 GeV). 

The total branching ratios for stop decay into the neutral 
channel, i, ; — > tx% 2 and tne charged channel, tj — > b\i, 
depend only on MSSM Yukawa couplings and phase space. 
Provided the mass splitting between stops and higgsinos is 
large enough that the top-higgsino mode is open, the stop 
branching ratios into charged and neutral higgsino modes are 
comparable in magnitude over the large majority of parameter 
space. This leads to a sizeable pair-wise branching ratio into 
the mixed final state, tt* — > tbx°x ± - Fig. [I] shows stop pair 
BRs to both the mixed channel and the usual search channel 
tt + 2x° in the mj-cos^t plane. The total erx BR into the 
tt + I$t fi nai state can easily be suppressed by a factor of two, 
while the erx BR into the mixed mode is often comparable 
or larger. The sbottom-like signal, tt* —> bb + 2x ± , makes 
up a distant third, except in regions where the neutral decay 
channel is suppressed by available phase-space. 

Assuming the gluino is inaccessible, the most stringent col- 
lider bounds on this minimal natural SUSY spectrum are the 
LEP chargino bounds |0. We assume for simplicity that 
the nondegenerate LEP bound m x ± > 103.5 GeV applies, 
though for sufficiently small splittings between the chargino 
and the (N)LSP the limits could be mildly relaxed. Limits 
from direct stop searches [4 5 1 have some reach in the region 
of interest (m x o > 100 GeV, > m x o + m t ) but still leave 
most of the region unconstrained, particularly when the signal 
is reduced by a branching ratio BR(tt — > tt + 1$t) ~ 0.5. In 
principle, the CMS razor analysis [6 | could also place inter- 
esting limits if analysed subject to this signal hypothesis. 

In what follows, we use a reference Higgsino sector 
where the mass splitting is taken to be (m x o , m x ± , m x o ) = 
[m x o , m x o + 6 GeV, m x o + 12 GeV). This splitting is large 
enough to reflect the possible influence of nearby gauginos but 
small enough that all higgsinos produced in the decay of 300- 
600 GeV stops appear primarily as missing energy: the addi- 
tional soft daughters produced in the decay of heavier higgsi- 
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FIG. 1: Stop pair branching fractions into (left) the mixed mode 



t + X°,b + x*\ (right) tt + 2 X ' J , for tan/3 
Higgsinos with mass /i = 100 GeV. 



20 and degenerate 



nos down to the (N)LSP x? are not boosted enough to pass se- 
lection cuts for hard isolated objects. We use a reference stop 
pair branching fraction of BR(tt* — > tb+tfr+X) = 0.5, and 
show results for a reference signal point with m t ~ = 500 GeV, 
m x o = 200 GeV, Xi = unless otherwise specified. In all 
signal points we take t\ — cos8 t iji + sin Q t l^ with mixing 
angle cos 6 t = 0.833. 

We have checked that the modes t — > t\\ and t — > tx\ — > 
txl + X pass our selection cuts with efficiencies differing by 
O(10%), with the single largest difference in efficiency com- 
ing from isolated lepton acceptance. Thus to good approxima- 
tion the asymmetric signal is insensitive to the details of the 
neutralino mixings. We have chosen a relatively large mass 
splitting for our reference Higgsino sector; in more degenerate 
spectra, the difference between x\ an d X2 w °uld be negligibly 
small. 



TOPNESS 

The dominant backgrounds to semileptonic stop searches 
are dileptonic ti where one of the leptons is either too soft or 
too forward to be identified, and tt events with one lepton and 
one unidentified r. Although three of the six partons in these 
final states are missing or unidentified, these backgrounds still 
contain a large amount of kinematic information which can be 
used to identify events consistent with top quark pair produc- 
tion. 

Much literature has been devoted to kinematic variables 
which can identify particle masses in the presence of multiple 
invisible particles. Two of the most studied variables are the 
stransverse mass, Mt2 0, and the contransverse mass Mqt 
[ 8 1 . Both of these variables admit straightforward extensions 
to the asymmetric decay chains that appear in top backgrounds 
with missing leptons, as was studied for Mt2 m 0- 

We propose here a novel alternative. Dileptonic top events 
are reconstmctible when both leptons are identified: the mass 
shell conditions provide enough constraints to completely 
solve for the unmeasured components of the neutrino mo- 
menta, up to discrete combinatoric and quadratic ambiguities. 



Once one of the leptons is lost, this is no longer true: the miss- 
ing particles are (by assumption) now a neutrino and a W, and 
one of the mass-shell conditions is lost along with the lepton, 
leaving an under-constrained system. 

We replace the missing mass-shell condition with the con- 
dition that the reconstructed center-of-mass energy of the 
event be minimized. As the PDFs fall off steeply with </s, this 
provides a good approximation to the true event kinematics. 
We construct a function S which quantifies how well an event 
can be reconstructed subject to the dileptonic top hypothesis: 
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where in the last term the sum runs over all 5 assumed fi- 
nal state particles. We have imposed transverse momentum 
conservation as well as the mass shell conditions p 2 u = 0, 
p w — rriyy to fix Ew, E v , p vx , and p vy in terms of the four 
remaining undetermined variables. The denominators at de- 
termine the relative weighting of the mass shell conditions, 
and should not be smaller than typical resolutions; we take 
aw = 5 GeV, at = 15 GeV, and acM = 1 TeV The value of 
S at its minimum quantifies how well an event can be recon- 
structed according to the dileptonic top pair hypothesis. The 
inputs to S are two jets, a lepton, and the f>T- To find the 
best possible reconstruction, we sum over both possible pair- 
ings of jets with reconstructed W bosons and keep the pairing 
which minimizes min S. When the event contains two identi- 
fied 6-jets, we use them as input to S; when the event contains 
only one identified b, we consider the two hardest untagged 
jets with |?7| < 2.5 and px > 20 GeV, and use the pair (b,j) 
which yields the minimum value for min S. We define top- 
ness as 



t = m(min S). 



(2) 



Minimization of S is a nontrivial computational problem. In 
our implementation we use 10 iterations of the Nelder-Mead 
algorithm per event. In general this is not sufficient to find 
the global minimum; however, it will find a minimum that 
is sufficiently close to the global minimum that cuts and dis- 
tributions are insensitive to any difference. We show distri- 
butions of topness for the major dileptonic and one-^-one-r 
top backgrounds as well as signal in Fig. [2] In the left panel 
of Fig. [3] we compare the performance of topness to both the 
asymmetric implementation of Mct which we find most ef- 
fective, and Mt2, the Mt2 variant identified as most effective 
in (9|. Events shown here have passed preselection cuts as 
described in the text below. 

Our present interest is in the asymmetric stop decay mode, 
but we emphasize that topness is useful in any search where 
the background is dominated by dileptonic tops with a missed 
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FIG. 2: Left: unit-normalized topness distributions for events passing 
preselection cuts as described in the text (signal, red, solid; dileptonic 
top, blue, dashed; one I, one r top, cyan, dotted). Right: topness dis- 
tributions for the dileptonic background broken down into samples 
with two truth b jets (blue, dashed) and one truth b jet (purple, dot- 
ted). 




FIG. 3: Performance comparison of topness (red, solid) to the vari- 
ables Mct (blue, dotted) and M^ 2 (purple, dashed) for the asym- 
metric signal (left) and the symmetric signal ti + 2\i (right). The 
quantity plotted is the gain in signal significance, assuming Gaussian 
statistics, as a function of signal efficiency. Results are shown for 
m; = 500 GeV and m^o — 200 GeV in both channels. 

lepton, most notably stop searches in other channels. The 
right panel of Fig.[3]shows the relative performance of topness 
on the traditional stop signal, tt* — > ti + 2xi, followed by 
semileptonic top decay, with backgrounds given by missed- 
lepton ti + 2j events. Events shown here have passed ob- 
ject selection cuts analogous to those in JT0|, including the 
requirement that mr(l,^T) > 150 GeV. Topness performs 
comparably to Mjf 2 m me symmetric search channel, realiz- 
ing (by a slim margin) the largest gain in significance among 
all three variables. The dilution in the efficacy of topness in 
the ti + 2\° channel, relative to the tb + B~t channel, is be- 
cause the extra jets in the events give more possibilities for 
signal to accidentally reconstruct as a top-like event. We do 
comment, however, that our implementation of topness was 
not optimized for the ti + I$t signal. 

HUNTING ASYMMETRIC STOPS 

In this section we discuss a search strategy for stop pair 
production in the mode tt — > tb+ MET. We target semilep- 
tonic top decay, so the final state of interest contains one 
£, 2 6-jets, and missing energy. We impose the following 
cuts at preselection: exactly one lepton satisfying p^ > 20 
GeV, p e > 25 GeV, \<q t \ < 2.5; rf T > 200 GeV, with 
mr(l,ftr) > 150 GeV; and at least two jets with pr > 20 
GeV and \t]\ < 2.5, at least one of which must be 6-tagged. 



The cut on wit suppresses all backgrounds where the J^t 
arises from a single W, in particular semileptonic tt and the 
enormous W+ jets, which is further suppressed by the 6-tag 
requirement. The major remaining background is therefore 
dileptonic top pair events where one of the leptons is not iden- 
tified, either because it falls outside acceptance, or because it 
is a non-identified r. A secondary background is the asso- 
ciated production of a top with a W boson, again with dou- 
bly dileptonic decays and a missed lepton. All major SM 
backgrounds can be reduced by identifying softer leptons in 
the event; we thus reject events containing identified hadronic 
taus with 1 77 1 < 2.5 andpx > 20 GeV or additional (isolated) 
leptons with \r]\ < 2.5 and pt > 15 GeV. Importantly, the 
additional soft decay products of the heavier Higgsinos in sig- 
nal events have negligible impact on the ability of signal to 
pass the veto. More aggressive vetos, as in O, would further 
reduce the backgrounds at little cost to signal. 

In addition to cuts on the hardness of final state particles 
and I$t, we add three novel variables which discriminate sig- 
nal and background. First and by far the most important is 
topness, discussed in the previous section. Another useful 
variable can be constructed by exploiting the asymmetric sig- 
nal kinematics. The 6-jet coming from the decay t — b\ + is 
typically much harder than the daughters of the top. The pt 
asymmetry 



PTb x + Pre 

of the lepton and the highest pr fo-jet, is thus useful for distin- 
guishing signal and background. We also employ a centrality 
variable, 

C = max (\Ar) juj2+e+fT | , \A>r)j2,ji+t+? T \) > ( 4 ) 

formed from the two highest px jets j% and j'2 in the event as 
well as the lepton and missing momentum. Centrality is typ- 
ically larger for backgrounds than for signal, reflecting both 
the larger signal masses and the different kinematics of top 
versus stop pair production ifTTIl . Distributions of r PT and C 
are shown in Fig. |4] This particular cut on rapidity separations 
is most useful when used in conjunction with topness, because 
background events with large topness have often selected an 
ISR jet in place of one of the b jets, and this ISR jet is distinct 
in rapidity from the other objects in the event. 

We normalize signal lfT2l and background ti [ 1 3 1 and 
tW 04 1 processes to inclusive NLO+N(N)LL 8 TeV cross- 
sections. Events are generated in Madgraph [15], showered 
in Pythia lfl6l . and clustered in FastJet using the anti-fc^ al- 
gorithm with R = 0.4 1171 . In generating tW + Ij events, 
we forbid ti events from contributing when the momentum 
in one of the internal top propagators lies in the window 
\p 2 — m-t\ < 15r t lfl8l IT9l . Leptons are declared isolated if 
the scalar sum-p^ deposited in a cone of radius Ri SO = 0.2 
around a lepton is less than r,; so = 0.2 times the lepton 
Pt- The isolation threshold is thus 4 GeV for a lepton with 
Pt = 20 GeV, which matches well with experimentally em- 
ployed criteria at threshold. We identify hadronic rs when 
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FIG. 4: Unit-normalized distributions of the variables r PT (left) and 
C (right), for events passing preselection cuts as described in the text. 
Signal is red and solid; dileptonic top is blue and dashed; one I, one 
r top is cyan and dotted. 
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TABLE I: Cut flow for the reference signal point, with = 500 
GeV, m x = 200 GeV, and assuming a 50% branching ratio into 
the tb + 1$t final state. All cross-sections are measured in fb, and 
significances are shown for £ = 20 fb -1 at 8 TeV. 



the scalar sum-py deposited in a cone of radius Ri SO = 0.2 
around the visible hadronic r is less than r iso = 0.6 times the 
r px- These criteria are tuned to yield an efficiency of approx- 
imately 50%, comparable to experimental tau working points, 
on hadronic ts with px > 20 coming from the W in ti events. 
Finally, we assume a uniform probability Vb = 0.7 to tag a b- 
jet with pt > 20 GeV and \r/\ < 2.5. We then generate more 
than enough events to ensure our final results are insensitive 
to statistical fluctuations in our Monte Carlo simulations. 

In Table U we show the detailed cut flow for our reference 
signal working point as further specialized cuts are imposed. 
As the total number of expected signal and background events 
is small, we use Poisson statistics to evaluate signal signif- 
icance. Table [II] shows optimal cuts for other signal mass 
points. We find in general that the region of best significance 
is relatively flat as a function of the cuts, that is, a broad range 
of cut combinations yield similar final significances. 
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TABLE II: Signal significances and best cuts for £ = 20 fb~ at 8 
TeV. All masses and energies are in GeV. For the signal point (m^ = 
400, m x o = 100), the quoted significance is Gaussian. 



To conclude, stops in natural SUSY spectra will generi- 
cally have large branching fractions to the mixed final state 
tt* — >• tb + I[it + soft, substantially diluting the rate into the 
traditional search channel ti — > ti + If?- Despite the large ti 
backgrounds in the tb + Ifx final state, good search reach is 
achievable with the use of the novel variables introduced here, 
recovering LHC sensitivity to the natural region of parameter 
space. The novel variable topness is critical to this search and 
is useful in more generic contexts, matching or outperforming 
such kinematic variables as M^f 2 and Mct- 
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